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Abstract 

A microscopic treatment of fundamental edge magnetoplasmons (EMPs) 
along the edge of a double quantum well (DQW) is presented for strong mag- 
netic fields, low temperatures, and total filling factor v = 2. It is valid for 
lateral confining potentials that Landau level (LL) flattening can be neglected. 
The cyclotron and Zeeman energies are assumed larger than the DQW en- 
ergy splitting V A 2 + AT 2 , where A is the splitting of the isolated wells and T 
the tunneling matrix element. Using a random-phase approximation (RPA), 
which includes local and nonlocal contributions to the current density, it is 
shown that for negligible tunnel coupling TT <C A the inter-well Coulomb 
coupling leads to two DQW fundamental EMPs which are strongly renormal- 
ized in comparison with the decoupled, single-well fundamental EMP. These 
DQW modes can be modified further upon varying the inter-well distance d, 
along the z axis, and/or the separation of the wells' edges Ay along the y 
axis. The charge profile of the fast and slow DQW mode varies, respectively, 
in an acoustic and optical manner along the y axis and is not smooth on the £q 



scale. For strong tunneling A < 2T these DQW modes are essentially modi- 
fied when A is changed by applying a transverse electric field to the DQW. 
PACS 73.20.Dx, 73.40.Hm 

I. INTRODUCTION. 

Recently the interest in EMPs has been increased considerably [|TJ- [ 171 ; for a discussion 



of older studies see, e.g., Ref. Initially EMPs were theoretically considered only within 
essentially classical models ||, |[9 111 which the edge does not vary but the charge density 
profile at the edge does. Except for the complications due to the long-range Coulomb 
interaction, this edge-wave mechanism is the magnetic analog of that for the Kelvin wave 
T8[ , i.e., the chiral wave that occurs in tides at the edge of a rotating "shallow" sea, the 



angular velocity of the sea being the analogue of the cyclotron frequency. More recently a 
distinctly different, fully quantum mechanical edge- wave mechanism for EMPs was provided 
by the models of Refs. ||- and [[□]], |12]] . In these models typically only the edge position 
of an incompressible two-dimensional electron gas (2DEG) varies; with respect to that the 
density profile is that of the undisturbed 2DEG ||16|| . 

In Refs. [[]1|- [17 we proposed a quasi-microscopic model that effectively incorpo- 



rates both edge-wave mechanisms mentioned above. In this model we treated EMPs 
oc A(lj, q x , y) exp[— i(ut — q x x)], for a single quantum well with v = 1 and very low temper- 
atures T, for which the unperturbed density profile drops sharply at the edges on a length 
of the order of magnetic length £ . Such a profile is valid for ksT -C hv g /£o, where v g is the 
group velocity of the edge states, and the unperturbed electron density n (y), normalized 
to the bulk value no, is calculated as n (y)/no = {1 + &[(y re — 2/)/A)]}/2, where y re is the 
coordinate of the right edge and Q(y) the probability integral. However, the current density 



j^u, q x ,y) was obtained [|T7| when the components of the electric field of the wave, E, are 
smooth on the £q scale. This is not well justified. In addition, nonlocal effects and the 
screening of the 2DEG were neglected. Later on, within a RPA framework we showed [TP 
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that the effect of nonlocality and screening by the 2DEG on the fundamental EMP is small 
and justifies well, for weak dissipation, the results of Ref. |I7J for the fundamental EMP. 
Here, using the same RPA framework we consider the fundamental EMPs at the edge of a 
DQW for v — 2, when two occupied Landau sublevels (LSLs), with n = and the same spin 
orientation, correspond to different DQW quantum numbers a z = — and a z = +. As all 
studies we are aware of treat EMPs in single quantum well systems, the present microscopic 
generalization of the model of Refs. [|16|] and Jl7| to EMPs in DQWs appears amply justified. 

In Sec. II we present the one-electron characteristics of a DQW channel, the wave 
charge density, inhomogeneous along y on a scale ~ £q <C A, and the corresponding electric 
potential. Further, we derive the integral equations for the fundamental EMPs in a DQW 
for negligible and strong tunnel coupling. In Sec. Ill we obtain the corresponding dispersion 
relations of the fundamental EMPs for very low temperatures. Finally, in Sec. IV we discuss 
our theory and make concluding remarks. 



II. BASIC RELATIONS 

A. DQW channel characteristics 

We consider a DQW, of width W and of length L x = L, in the presence of a strong 
magnetic field B along the z axis. In the absence of a self-consistent wave potential V(r, t) = 
V(u , q x , y, z) exp[— i(u Q t — q x x)} + c.c, we assume that the motion of an electron in the 2D- 
plane (x,y) is decoupled from it's motion along the z axis. The latter is described by the 
Schrodinger equation 

h z X az (z) = h^-^ + V z ]X az (z) = E az X az (z)- (1) 

The DQW consists of two square wells along the z axis, with finite depth U, widths di and d r , 
and separated by a barrier of width d. In addition, we assume that m* is independent of z. 
As is usual in DQW studies, we assume that only a pair of tunnel-connected closely located 
levels, formed from the ground states of the left (I) and right (r) isolated QWs, is essential. 
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For small barrier penetration |2(J we can use these isolated r and I well eigenfunctions, 



denoted, respectively, by X r (z) and Xi(z) as a basis for solving Eq. ([]]). With a z = ± the 
corresponding normalized orbitals X az read p( 



X+ = N[Xi + (2T/(A + A T ))X r ], X_ = iV[X r - (2T/(A + A T ))X,]. (2) 



Here N = \J(A + At)/2At is the normalization factor, T = — Toexp(— nd) (TJ 
2ii 2 h 2 / m* nd^ 2 d^ 2 ) is the tunneling matrix element, and /t _1 ~ fi/\/2m*U . Further, A 
is the energy splitting of the isolated I and r levels, and A^ = V A 2 + 4T 2 is the splitting 
between the a 2 = + and a z = — levels; A^ takes into account tunnel branch-off. The 
eigenvalues E + and E- are +At/2 and —At/2, respectively. Here weak barrier penetration 
means exp(—2nd) <C 1 and nd^ r ^> tt. In the calculations of pertinent matrix elements we 



can approximate X^ r (z) by y2/di jr cos(ir(z — z^ r )/d^ r ) for \z — zi >r \ < di tr /2, and by for 
\z — Zi >r \ > di >r /2, where zi, and z r are the centers of I and r QWs, respectively. That is, 
the matrix elements can be calculated without taking into account the barrier penetration 
of the isolated well orbitals, cf. Ref. [p0[ . 

The lateral confining potential is taken flat in the interior of the 2DEG (V y = 0) and 
parabolic at its edges, e.g., V y = m*Q 2 (y — y r ) 2 /2, y > y r . V y is assumed smooth on the scale 
of £o = {h/m*uj c ) l l 2 such that Vt cu c , where tu c = \e\B/m* is the cyclotron frequency. We 
will also consider the case of a tunnel-decoupled DQW, with different confining potentials in 
the I and r QWs, i.e., = m*tt 2 (y - y i l l ' r) ) 2 /2, y > y^ with y<?) ^ y^. In the Landau 
gauge for the vector potential A = {—By, 0, 0) and without the interaction V(r, t), the total 
one-electron Hamiltonian hP is given as h° = h 2 d + h z , and 

h 2d = [{p x + eBy) 2 + f y ]/2m* + V y + g*fi B S z B/2 . (3) 

Here p is the momentum operator, g*(< 0) the effective Lande g-factor, and \xb the Bohr 
magneton. S z is the z-component of the spin operator with eigenvalues a — 1 and a = — 1 for 
spin up (|) and down (!), respectively. The eigenvalues and eigenfunctions corresponding 
to Eq. (|3|) for the right edge of the channel, where y = yo(k x ) = t\k x > y r , are well 
approximated |2l| by 



E a2d = E n>(T {k x ) = (n + l/2)hu c + m*Q 2 (y - y r ) 2 /2 + ag*p: B B/2, (4) 

and 

respectively. Here p = {x, y}, k r = y r /^o, «2d = {n, k x ,a}, ^ n (v) is a harmonic oscillator 
function, \a >= = 5 aai is the spin-wave function, and o\ = ±1. Then the eigenvalues 

of h° are E a = E a2d + E Qz and its eigenfunctions \a >= ip a2d (p)X az (z)\<r >. 

We will consider only the case of v — 2, when only two tunnel-split levels are occupied. 
These Landau sublevels (LSLs) have the same Landau and spin quantum numbers, n = 
and cr=l, but different quantum numbers a z = ±. In line with Refs. |^2| and [^3], we 
assume that hu c > At and \g*\p B B > At- Then the one-electron energy spectrum (for 
Uo > Ur) for the two occupied LSLs is given by 

E Q , kx ,i,± ee £ ,±(fc*) = [^c + m*ft 2 (?/o - - \9*\»bB ± A T ]/2. (6) 
Equation (^) leads to the group velocity of the edge states vf = ^ -1 9^o,±(^r + k^)/dk x = 



UVPk®^ /m*oj 2 with characteristic wave vector k®^ = (uj c /hfl)y2m*Ap Q , where A^ = 
Ep — [hu c — \g*\p B B ± At]/ 2 > 0, and E F is the Fermi energy. The edge of the LSL is 
denoted by yf = y r + l^k®^ = l^k^, where kf Q = k r + fc° ,:t , and = 2^. We can also 



write f s0 = E e0 /B, where E e0 = Vt^j2m* Ap /\e\ is the electric field associated with the 
slope of the confining potential V y at yf Q . 

For definiteness, we take the background dielectric constant e to be spatially homoge- 
neous. Assuming I^K^ 3> 1, we can consider an EMP along the right edge of the channel, 
of the form A(u, q x , y, z) exp[— i(ut — q x x)], totally independent of the left edge. 

B. Wave charge density and electric potential at the DQW edge 
As in Refs. [IDT] , [F2TJ , we assume that without interaction the one-electron density matrix 



P 



(o) is 

diagonal, i.e., < a\p^\(3 >= f a 5 a p, where f a = 1/[1 + exp((E a — EF)/k B T)] is the 



Fermi-Dirac function. In applying the RPA we follow the self-consistent field approach which 
is discussed, e.g., in Ref. [p^j . The one-electron Hamiltonian is H(t) = h° + V(x,y,z,t) in 
the presence of the interaction taken as a self-consistent wave potential V(r,t). Then the 
equation of motion for the density matrix p reads 

in^ = [H(t),p)-^{p-p^), (7) 

where [, ] denotes the commutator. On the right hand side (RHS) of Eq. (|7|) we have 
introduced phenomenologically the infinitesimal term oc 1/r (r — > oo) that leads to correct 
rules for passing an integration contour around the pole singularities, etc., cf. Refs. |21| - 



Z5fl . Notice that r — > oo corresponds to the collisionless case whereas a finite r provides 
the possibility of estimating very roughly the influence of collisions. The most effective 
scattering is related to intra-LSL and intra-edge transitions [17], PT . 



Applying the Laplace transformation over time to Eq. (0) and writing R(u) = e lujt pdt 
and R a p{uj) =< a\R(u)\/3 > we can present the solution of Eq. (|7|) as a power series in V 

oo 

R af} (u) = £ R^\u), (8) 

m=0 

where R^l(u) = (ifa/w)5ap- Because we consider linear EMP's, it is sufficient to take into 

account only the first two terms m = and m = 1 on the RHS of Eq. (H). Then in 

V(x, y, z, t) we can consider only the term V(u , q x , y, z) exp[— i(u t — q x x)], which leads to 

\ _ Kfp ~ fa) < u\V(u ,q x ,y,z)e iq * x \p > 
a(3[UJ) (LU-u )[E^-E a + nu + i?i/T} ' 1 ' 

Because V(r, t) conserves spin the spin index in Eq. is the same, a a — ap — 1, and 
therefore it can be dropped. For simplicity we will use the notation a = {n, k xa , a z }. Then 
taking the trace of the density matrix p with the electron charge density operator, eS(r — f ), 
we obtain the wave charge density as 

6p(t,x,y,z)=p(t,x,y,z) = — / dwe^^^SM^W^W, ( 10 ) 

Z7T J-oo+ir/ a p 

where rj > and ip a { r ) = i>a 2 d(p)Xa z {z)- From Eqs. (|) and ( p~0|) it follows that p(t, x, y, z) = 
p(t,q x ,y, z) exp(iq x x). Moreover, for t/r ^> 1, when contributions related to transitional 
processes are negligible, Eqs. @ and (^0|) lead to p(t, q x , y, z) = p(ui , q x , y, z) exp(— iui t). 



To simplify the calculations we assume that cZ; — ^ and d r — > 0, i.e., we neglect the 
thickness of both QWs, cf. Ref. |22|. This leads to X 2 r (z) 8{z ± d/2) for d lr -> 0. 



Then Poisson's equation gives the wave electric potential 4>(t, q x , y, z) induced by the density 
p{t,q x ,y,z) =p®(t,q x ,y)5(z + d/2)+pW(t,q x ,y)5(z-d/2) in the Z QW (z = -d/2) and r 
QW (z = d/2) as 

2 r°° 

<p(t,q x ,y,z = -d/2) = <f>i(t,q x ,y) = - dy'K lr (y,y';t } q x ) } (11) 

€ J —oo 

and 

2 f°° 

<f>(t, q x , y, z = d/2) = <f) r (t, q x , y) = - dy'K H {y, y'; t, q x ), (12) 

6 J —oo 

where 

K lr (y, y'; t,q x ) = K {\q x \\y - y'\)p {l) {t, q x , y') + K {\q x \d w )p^{t, q x , y'). (13) 
Here K (x) is the modified Bessel function and d yy i = \J(y — y') 2 + d 2 . 



C. Integral equations for EMP's 

For t/r ^> 1, we can change p^\t,q x ,y') to p^(cuo,q x ,y') (j = l,r) on the RHS of 
Eqs. ( p"TD and (0). Then, omitting the common factor exp(— iu>ot), we should change 
(fii,r(t, q x , y) to 0; jr .(co>o, q x , y)- in the absence of an external potential we have V(ujq, q x , y, z) = 
e(j)(uJo,q x ,y, z). Then using Eqs. (|ITi|)-(|i"2]) and integrating Eq. (PH|) over z once in the 
neighborhood of z = —d/2 and once in that of z = d/2, we obtain a system of two integral 
equations for p^ l '(u;o,q x ,y) and p( r '(uo, q Xl y)\ if we omit the subscript in u they read 

p U) y) = E EE i //3 w^+^l z J*L n w fc ^> 

eL n a ,n p =0 k xa a zA Ep-E a + TlUJ + III / T 

/oo roo 
dy / dy' U na n (y, k xa , k x p){[< a z \(3 z > z K lr (y,y';u,q x )] + [r <-+ I]}; (14) 
-co J — oo 

where n nan/3 (y, fe^, fe^) = ^ nct {y ~ yo(k xa ))^ n0 (y - yo{k x p)), k x p = k xa - q x , < a z \p z >i= 

J-d/2-S dzX *z X Pz and <<*z\Pz>r= Id/2-0 dzX az Xp z . 
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The integral equations (|14]) take into account the nonlocality in the electron current 
density oc J dy'an V (y,y')E u (y'), and the screening by the edge and bulk states of the 2DEG. 
In Ref. |26| it is shown that screening by the edge states can be strong. We consider very 
low temperatures T satisfying hv g0 ^> loksT and use the long- wave length limit q x £o <C 1. 
Further we will obtain only the low-frequency solutions of Eqs. (Ill) for i) negligibly small 



tunneling w<w c and ii) strong tunneling uo -C Ay/ft. In both cases we assume u c > At/%. 

Case ii). Comparing the terms oc fp* on the RHS of Eq. (|14[) , for given rip* = and 
(3* = + or — , we find that the contribution to the sum over n a and a z with n a = rip* = 
and a z = /?*, is much larger than that of any other term in this sum (e.g., the term with 
n a = and a z ^ (3*) or that of all other terms with n a > 0. The small parameter is 
h\u - q x vf (k x p)\/A T < 1, where vf (k x/3 ) = h^dE^k^) / 'dk x(3 is the group velocity of 
any occupied state {0, k x p, of the f3* z LSL. We assume that both occupied LSLs have one 
intersection with the Fermi level at the edge of the channel and denote the group velocity of 
its edge states by vf = vf (k^ ). The small parameter given above implies hq x Vg /A T <C 1, 
since v ~ is typically larger than v£ Q . Similar results hold for the terms oc f a * in the sum over 
rip and (3 Z . Hence for Tiuj <C A T and hq x vf <C A T the terms with n a ^ rip or/and a z ^ (3 Z 
can be neglected. This leads to a system of simpler integral equations for p")(u>,q x ,y) and 
p^(u,q x ,y) 

Jj)(,. n v ) — ^_ V - a ^ Qz >j (f°> k ™-i*> a z ~ fo^ xa ,a z ) tt f ?/ u u _ n \ 
eL E az (k xa - q x ) - E 0)Otz (k xa ) + huj + %h r 



x 



f dy f dy' U 00 (y,k xaj k xa - q x ){[< a z \a z >i K lr (y,y';uj,q x )} + [r ^ /]}. (15) 



Case i). For negligibly small tunneling we have 2|T|/A -C 1 and from Eq. (0) it follows 
that X + (z) m Xi(z) and X-(z) ~ X r (z). Again we can take n a = rip = 0. Then Eq. ( HM 



for lo < lo c , q x v% < u c , j = I, and a z = j3 z = + gives 



p®(u,q a ,y) = — [°° dk xa — — /0 ' fc 7' fe F + ,{°' fc r; + , ^ . %l n 0Q (y,k xa ,k xa ~q x ) 

ne J-oo E 0:+ (k xa - q x ) - E 0+ (k xa ) + huj + ih/T 
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/oo roc 
dy / dy' Tl 00 (y, k xai k xa - q x ) K lr {y, y'\ u, q x ) ; (16) 
-oo J —oo 

for j = r and a z = (3 Z = — the result is 

pV(u, q x , y) = — f°° dk xa ^f^^fi — ^ n 00 (y, k xa - q x ) 

ne J-oo E ,-{k xa - q x ) - E „(k xa ) + hu + m/r 

/CO /"CO 
/ dy' n 00 (y, - iT ri (y, y'\ to, q x ). (17) 

-oo J —oo 

We will study only the fundamental EMPs following from Eqs. ([15[) - (|TTD . 



III. FUNDAMENTAL EMPS, v = 2 



A. Negligibly small tunnel coupling 



We now consider case (i), when 2|T|/A <C 1 holds, and we have the system of Eqs. 
( P^D and ( [T7| ) for p®(u, q x , y) and p^(uj, q x , y) coupled only by the Coulomb interaction. In 
the assumed long-wavelength limit q x £o < 1 we can approximate fo,k xa -q x ,± ~ fo,k xa ,± by 
QxS(k xa — kf ) and neglect the small shift oc q x £ in the argument of ^o(x — yo{k xa ) + l\q x )-, 
i.e., we can approximate the latter by *$?o(x — yo(k xa )). Then from Eqs. fll6[) and CflTf) , after 
integration over k xa , we obtain 



P il) (uj,q x ,y) = 4" g '* g(y + f°! I" dy r dy' ^{y - y+) ;/:,:.,,. ) (IS) 

7me uj — q x vJ Q + i/t J-oo 



oo roo 



—oo 



and 



OO ("OO 



P (r) (^fe2/) = 4 VrQ ] I" d y I" dy' <{y-y7v) K rl (y,y';u, qx ), (19) 



nhe lu — q x v g0 + i/r j -oo j-oo 
respectively. The solution of Eqs. ([18]) and ([19]) can be written as 

p«\u, q x , y) = p®(u, q x ) *t(y - y+), p^(co, q x , y) = p<r\u, q x ) * 2 (y - y ; ). (20) 

Since X+(z) ~ Xi(z) and X_(z) ~ X r (^), we can consider v£ Q and y^ as parameters of the 
/ QW that are independent of those of the r QW v~ (y~ ) ; the same holds for E 0)+ (k x ) and 
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Eo,-(k x ). Thus, in Eqs. ([1S|)-(|20D we can take, e.g., Ay = y~ — y^o ^ ^o- I n particular, here 
we can assume different confining potentials in the I and r QWs, Vj l ' r \ 

Substituting Eq. (p0|) into Eqs. (|18|),(|i~9|) gives a system of two linear homogeneous 
equations for the unknowns p^ l ' r \uj, q x ). The condition of a nontrivial solution gives the 
dispersion relation (DR) for the renormalized, by the Coulomb coupling, single-well funda- 
mental EMP (w®(q x ) = wW) as 



2 2 

"2 = q*V 9+ -- + q x a 00 (q x ) —r- - a z q x [(^-) 2 b 2 0Q (q x , Ay, d) + V 2 _] l '\ (21) 
where V g ± = (v~ ± v g0 )/2 and a 2 = ± corresponds to the ± LSL. Here 

/oo roc , 

dx / dx' *l{x)*l{x') K Q {\q x \J[{x-x') + Ay\* + d?) (22) 
-oo J —oo 

and aoofe) = b 00 (q x , 0,0) « [h^l/l^o) + 3/4] [T|; it is assumed that v g0 > v+ > 0. 
Without inter-QW Coulomb coupling, i.e., with boo(q x , Ay, d) = in Eq. fl21"|), we obtain the 
DR of the fundamental EMPs for totally decoupled QWs as cj'l ^ = q x v~ + [e 2 / irhe]q x aoQ(q x ) — 
i/r, for the r QW, and as lu+ = q x v g0 + [e 2 / TThe\q x a 00 (q x ) — i/r, for the I QW. If we neglect 
dissipation, this result coincides with that of Ref. |T7| for v — 1 in a single QW. For a 
GaAs-based DQW (e = 12.5) we have 2e 2 /irTie ~ 10 7 cm/sec which is essentially larger than 
the typical v g0 < 10 6 cm/sec , cf. [|7|], ]17|. Moreover, for D = (Ay 2 + d 2 ) 1 / 2 < 0A/q x > £ 



it follows that b 0Q (q x , Ay, d) > 1. Then for q x D < 1 and D 2 /£ 2 > 1 Eq. fl2l]) gives the DR 
of the renormalized fundamental EMP of the r QW, i.e., the fast EMP, as 



(<0 



^ = q x V g+ + q x A~ [2H-TTT) ~ Ht) + °- 85 ] " l \ ( 23 ) 
and that of the / QW, i.e., the slow EMP, as 

4 C) = QxV g+ + M^) + 0.65] - i. (24) 

Substituting cj_ and given by Eqs. (53) and (|2H|), into Eq. (JTM) we obtain 



p {l \uj-\q x )/p {r \u-\q x ) ~ 1 andpW(wP,g a .)/pW(i<;P,g a .) « -1, respectively. Forv 90 = v+ 
these ratios are exactly 1 and —1, respectively. Then, as can be seen from Eq. (P0f), the 
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renormalized fundamental EMP of the r QW, Eq. ( p3|) , shows an acoustic spatial behavior 
in the y direction and that of the I QW, Eq. (p4]), an optical one though its DR is purely 
acoustic. 

For negligible dissipation the DRs ( p3[ ) and (|24]) are shown in Fig. 1 by the solid [D 2 j£\ = 
3) and dashed [D 2 j£\ = 10) curves, respectively, for t>~ rs = 10 6 cm/sec, e = 12.5, and 
uj* = e 2 /nhe£o. The dotted curve shows the DR of an isolated, single-well fundamental EMP 
ur_ = oj+ = q x Vg + [e 2 /Trhe}q x [ln(l / (\q x \£ )) + 3/4]. The upper solid and dashed curves 
represent the fast EMP and the lower ones the slow EMP. 



B. Strong tunnel coupling 

We now consider case (ii) in which 2\T\ > A holds. As in Sec. Ill A, we make the 
approximations f ,k xa - qx ,± ~ fo,k xa ,± ~ Qx5(k xa - kf Q ) and ^I (x - VoiKa) + ^Qx) ~ ^o(^ - 
yo(k xa )). Then integration over k xa in Eq. ( p3|) results in 

/oo roc 
dy / dy' y 2 (y - y%){[< a z \a z K lr (y,y';uj,q x )] + [r «-> I]}, (25) 
-oo J — oo 

where j = l,r. The solution of Eqs. (p5) can be written as 

p®(w,q x ,y) = ]T < a Mz >i Pa.(w,qx) *o(</ -<o), (26) 

a z =± 



where i — l,r. Substituting Eq. ( f26| ) into Eq. (p5| ) and demanding that the factors in front 
of ^l(y — yto) on both sides of Eq. (^) be equal for j = I and for j = r, leads to a system 
of two linear equations for p±(u, q x ) written as 

P«M qx){^-(u - q x v% + z/r) - iV 4 [(l + M- 8 )a o(&) + 2Ar 4 & o(&, 0, d)]} 

^ q x 

-p_ az (LU,q x )N 4 {[l + M- 8 }b 00 (q x ,Ay,d) + 2M~% (q x ,Ay,0)} = 0. (27) 
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Here N 2 = (A + A T ) /2A T , M 2 = (A + A T ) /2T, a* = - for the first equation and a 2 = + 
for the second one. The condition for a nontrivial solution of these equations gives the DR 
of the renormalized, due to tunneling and Coulomb coupling, single-well fundamental EMP 

,2 



(tc) 



q x V g+ -- + q x -4-iV 4 {[l + M- S }a 00 (q x ) + 2M~% m {q x , 0, d)\ 



e 



- a M^/ N8 ^ + M ' 8 ) b oo(^ A ^ d ) + 2M-% (q x , Ay, 0)] 2 + V g 2 _} 



1/2 



(28) 



Notice that for finite A and T — > 0, i.e., without tunnel coupling, Eq. fl28|) coincides with 
Eq. ( fnf ). Further, the first term under the square root in Eq. (p8|) is typically much larger 
than the second one. Then 



LO 



(tc) 



q x v g+ h g x 



aoo(q x ) - a z b 0Q (q x , Ay, d) 



2e 2 T 2 



7rhe A 2 ^ 



a>oo(<lx) ~ b 00 (q x ,0,d) 



(29) 



on the RHS we have neglected the cotribution a z [b 00 (q x , Ay, 0)—b 00 (q x , Ay, d)}, which is small 
in comparison with [a 00 (q x ) — b 00 (q x , 0, d)} due to the assumed condition Ay 2 /d 2 3> 1. Now 
in the actual region 1 > d/£o > 1/4 we numerically evaluate [aoo(q x ) — b o(q x , 0, d)} « d/£o 
with relative precision < 10%. Then Eq. (]29l ) gives the DR of the fast EMP as 

- 2 1 , . ,£> _ 2d T 2 , i 



LO 



(tc) 



q x V g+ + q x 



[21n(- 



and that of the slow EMP as 



lo 



(tc) 



ln( 4 ) + °- 85 -4 A 2 + 4T- r 



T 1 



(30) 



e 2 D 2d 



(31) 



where Ay/£ Q = £ [k° e - - k° e >+] 
and Apo 3> At/ 2. Using Eqs. 



c /fi)[(A 2 + 4T 2 )/2^ C A F0 ] 1/2 for A± = A F0 =f A t /2 



(gl]), and (^) we obtain p_(u;i* c) , q x ) / p+(u K l c> , q 



Xtc) 



1 



for the fast mode and p-(u+ , q x )/p + (u+, q x ) = — 1 for the slow mode. The corresponding 
charge density profiles are obtained from Eq. (EBf) as 
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P(0( "%; g '\ y) = N 2 {M-^l(y - y; ) ± *g(y - 
p-(u± \q x ) 



N 2 M(y ~ y; ) ± M-^l(y - y+)}, (32) 



where the upper (lower) sign corresponds to the fast (slow) mode. 

For strong tunnel coupling A < 2|T| = lmeV we plot in Fig. 2 the group velocities 
of both DQW fundamental EMPs, obtained from Eqs. ( j30|) and ([Hp, as function of A 
with = e 2 /nhe and A* = y^irieV. A can be changed experimentally by a transverse 
electric field applied to the DQW structure. Notice that the group velocity of the fast mode, 
du^/dq x , depends on q x whereas that of the slow mode, duj+ c \q x ) / dq X} does not. For a 
GaAs-based DQW, SI = 7.8 x 10 11 sec" 1 ^7|1, and B = 9 Tesla we have uJVt « 30, and 



Ai?o = 4meV. The value of d/£o is 0.5 for the solid and dotted curves and 1 for the dashed 
curves. The upper solid and dashed curves represent the group velocity of the fast EMP for 
q x £o = 10~ 3 and the lower ones that of the slow EMP. The dotted curve represents the group 
velocity of the fast DQW fundamental EMP for q x £ = 10~ 2 . For A = we have Ay/£ w 2.7 
(A T /2A F0 = 1/8 < 1) and Ay/£ w 5.4 for A = meV (A T /2A F0 = 1/4 < 1). Thus 
the assumed conditions (Ay/d) 2 ^> 1, etc., are well satisfied in Fig. 2. In addition, we have 
vf Q ~ 5 x 10 5 [1 =F A T /(4A^ )] cm/sec and the term q x (v~ + v^ )/2 is independent of A or 
T. As is clear from Fig. 2, the change of A/A* from to 1 changes the group velocities of 
both DQW fundamental EMPs appreciably and especially that of the slow one which shows 
a 1.78 times increase of its group velocity within that region. Finally, for larger d/£o the 
group velocity of the slow mode shows a stronger dependence on A (for small A) whereas 
that of the fast mode shows a weaker one. 

IV. DISCUSSION AND CONCLUDING REMARKS 

The above treatment has shown that in a typical DQW the Coulomb and tunnel coupling 
lead to two new DQW fundamental EMPs whose dispersion laws and spatial structure, along 
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the y-axis, are essentially different than those of single- well fundamental EMPs. These single- 
well EMPs follow from the former if we neglect inter-well coupling. The fast EMP behaves 
as uj oc q x \n.(l/q x ) whereas the slow EMP behaves as uj oc q x . The latter has an acoustic 
dispersion law though its spatial structure, along the y axis, is optical. Indeed, for both 
negligible and strong tunneling, when the DQW levels are in resonance and A vanishes, the 
charge localized at the edge y~ of the - LSL has exactly the opposite value of that localized 
at the edge y^ of the + LSL where Ay = y~ — y^~ ^> £ . The group velocity of the slow 
EMP is much smaller than that of the fast EMP. Moreover, they can be easily modulated 
by applying a transverse electric field to the DQW structure in the case of strong tunnel 
coupling. This renormalization of single-well fundamental EMPs is essentially different than 
that of the right- and left-edge fundamental EMPs in a single-well 2DEG channel with 
relatively short width, q x W -C 1. As shown in Ref. p] , in the latter case the renormalized 



EMPs have dispersion laws uj oc ±g x Wln(l / q x ) that give the same absolute value for a group 
velocity. 

We have neglected the possible spatial inhomogeneity of the background dielectric con- 
stant e along the z direction. In GaAlAs/GaAs QWs such an inhomogeneity is relatively 
small. Within the RPA we can include the effect of scattering only in a phenomenological 
manner through the small parameter 1/r. This effect can be treated approximately along 



the lines of Refs. |17| and ||21|| . We have also neglected the possible screening influence of 
the gates or other free charges outside the spacer layer. Their influence on the slow mode, 
if they are at a distance much greater than Ay far from y^r , should be negligible. We have 
also assumed that the confining potential, flat in the interior of the channel, is smooth on 
the magnetic length scale but sufficiently steep at the edges that Landau-level flattening 



|28|1 can be neglected [29. In addition, though the thicknesses d[ and d r of the QWs have 



been neglected, we believe that this often used approximation |22j does not affect our results 
essentially. 

Note added in proof. Remarkably the results of Sec. Ill A partly overlap with the 
M = 1 results of Ref. 30, obtained in the time-dependent Hartree-Fock approximation for 
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the particular M = 1 edge excitations of a double quantum-dot system with zero tunneling 
between the dots and total v = 2 in the limit of large dots. In this limit the dot radius R is 
large compared to both the magnetic length and the dot separation d. Then the energy 
of the in-phase collective magnetoplasmon excitations E^ pl is given by Eq. (29) of Ref. 30 
and that of the out-of-phase excitations E~ l by Eq. (30) of Ref. 30. The former follows 
from Eq. (23) and the latter from Eq. (24) of Sec. Ill A if we make the changes q x — > l/R, 
e — > 1 , and neglect q x V g+ and the important terms 0.85 in Eq. (23) and 0.65 in Eq. (24) . 
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FIGURES 

FIG. 1. Dispersion relation of the DQW fundamental EMPs, for negligible tunnel coupling 
2|T|/A <C 1. The solid curves are for D 2 jl% = 3 and dashed ones for D 2 jl\ = 10; uj* = e 2 / (trheio) . 
The dotted curve is the single-well fundamental EMP without inter-well Coulomb coupling. The 
upper solid and dashed curves represent the fast EMP and the lower ones the slow EMP. 

FIG. 2. Group velocity of the DQW fundamental EMPs for strong tunnel coupling, 
A < 2\T\ = lmeV, as a function of A; v* = e 2 /(irhe) and A* = \/3meV. The value of d/£o 
is 0.5 for the solid and dotted curves and 1 for the dashed curves. The upper solid and dashed 
curves represent the fast EMP, for q x £o = 10 -3 , and the lower ones the slow EMP, whose velocity 
is independent of q x . The dotted curve represents the fast EMP for q x £o = 10~ 2 . 
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